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Abstract 

By using a matrix technique, which allows to identify directly the ladder operators, 
the Penning trap coherent states are derived as eigenstates of the appropriate annihilation 
operators. These states are compared with the ones obtained through the displacement 
operator. The associated wave functions and mean values for some relevant operators in 
these states are also evaluated. It turns out that the Penning trap coherent states minimize 
the Heisenberg uncertainty relation. 

PACS: 03.65.Ge, 03.65.Sq, 37.10.Ty, 37.30.+i 

1 Introduction 

The coherent states (CS) approach to quantum physical systems [1-3] constitutes nowadays an 
alternative to the standard method, which address the same problem in terms of energy eigen- 
states and eigenvalues. Along the years the CS have been derived for plenty of Hamiltonians 
having either a ground or a top state, and some of them admit a group theoretical construc- 
tion in which this state is acted on by an appropriate displacement operator [2,4]. However, 
there exist interesting physical systems for which the Hamiltonians have neither ground nor 
top state [5,6], but it is required anyway a systematic technique to build up the corresponding 
CS. One of those systems consists of a charged particle in an ideal Penning trap [7, 8]. Such 
an arrangement, sometimes called Geonium atom, has been largely used to perform high pre- 
cision measurements of fundamental properties of particles [7]. Moreover, it could be used 
to test and/or control some intrinsically quantum phenomena as entanglement, decoherence, 
wavepacket reduction, etc [8-10]. 

In this paper we are going to address, from a coherent states viewpoint, the quantum motion 
of a charged particle in a Penning trap. With this aim, in section 2 we will present some gen- 
eralities of the standard coherent states. In sections 3 and 4 we will introduce the Penning trap 
Hamiltonian and discuss its corresponding algebraic structure. It will be shown that the system 
possesses a certain "extremal" state, which plays the role of a ground state although there is not 
a minimum energy eigenvalue. In section 5 we will construct the wavefunction associated to 
the extremal state, while in section 6 we will perform the corresponding CS construction. The 
mean values of some physical quantities in the CS will be calculated in section 7. Finally, in 
section 8 our conclusions will be presented. 
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2 Standard coherent states 



Glauber definitions of CS are based on properties of the harmonic oscillator [11], which have 
been applied to several different systems (see e.g. [1-3]): 

(1) The CS 1^) are eigenstates of the annihilation operator a: 

a\z) = z\z), z E C (1) 

(2) They arise from acting the displacement operator on the ground state, 

\z) = D{z)\i)o), D{z) =exp{za^ - z*a), (2) 

being the creation operator. 

(3) The CS satisfy the minimum Heisenberg uncertainty relation for X and P, 

(AX),(AP), = h/2, (3) 

where {AO)l = {z\{0 - = (C^)^ - {0)l is the mean square deviation of an 

observable O in the state l^;). 

It is worth to notice an additional property of the standard CS, which is relevant since some 
authors consider it as the fourth CS definition. It is the completeness relationship ^ J \ z){z\d?z = 
1, where 1 is the identity operator. In fact, the CS are overcomplete in the sense that for any con- 
vergent sequence of complex numbers z„ the corresponding CS \zn) form a complete set [12]. 

For systems different from the harmonic oscillator, these definitions lead to different sets of 
CS. In this paper we will use the first and second definitions to find the CS for a charged particle 
in an ideal Penning trap; we will show that they satisfy as well equation ([3]). 



3 Penning trap Hamiltonian and the matrix A 

Let us consider a spinless particle of mass m and electric charge e inside of an ideal Penning 
trap, i.e., under the influence of a constant homogeneous magnetic field pointing along the z- 
direction B = Bk, and a static electric field E = — V$(r), both arising from the following 
vector and quadrupole scalar potentials: 

A(r) = -l-fx B, $(f) = ^{x^ + y^- 2z^). (4) 

Throughout this paper, the small letters r, p, x, y, z, p^, Py, Px will denote either classical coor- 
dinates and momenta or the eigenvalues associated to the corresponding quantum operators, 
the last ones being represented by capital letters R, P, X, Y, Z, P^, Py, Pz- The Hamiltonian 
describing our quantum system is given by 

H = —{P--AiR)Y + e^R) = — + bLz + -[ib' + v){X' + Y')-2vZ^], (5) 

where L = R x P is the angular momentum operator, b = — v = and we take by 
simplicity 6 > 0. To ensure that the particle is trapped inside the cavity, some restrictions on the 
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parameters b, v have to be taken: first of aWv < in order that the 2; -motion is bounded (so that 
this mode is characterized by a standard oscillator Hamiltonian). However, the corresponding 
repulsive oscillators in the x — y plane do not have to destroy the trapped motion induced by 
the magnetic field, which is achieved by taking b"^ + v > 0. 

From now on we will assume that m — 1 and h — 1. Note that this assumption is equiv- 
alent to the following procedure: (i) first making the operator changes Ri — Risjm/h, Pi — 
P-J Vmh, i = 1,2,3, H — H/ h; (ii) then dropping the hats in order to simplify the notation. 
Thus, the Hamiltonian we are dealing with reads 



H = ^ + bL, + ^[{b' + v){X' + Y')-2vZ'], 



(6) 



iS 



where [Ri, Pj 

It is useful to work in the Heisenberg picture in which the evolution of the operator vector 
r]{t) = W{t)riU (t) is simply determined from a matrix equation: 



dr]{t) 
dt 



U\t)[iH,ri\U{t) = U\t)kr]U{t)^Kr]{t) 



(7) 



where r] = {R, P)^ involves the observables R, P in the Schrodinger picture, the superindex 
denotes to transpose the involved vector, U (t) is the evolution operator such that U (0) = 1. 
The calculation of [iH, 77] = At] leads to 



(8) 



Let us find next the right (u) and left (/) eigenvectors of the matrix A, which are called eigen- 
vectors and eigenforms respectively. Since A is non Hermitian, the eigenforms / are not neces- 
sarily the adjoints of the eigenvectors u. In order to determine both, we solve in the first place 
the characteristic equation of A: 
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|A - Al| - A*5 + - v{Sb^ + 3v)X'^ - 2v^ = 0. 
Thus, the eigenvalues are ±Ai = ±iuji, ±A2 = ±iuJ2: ±-^3 = ±ius, where 
= 6 + + V, U2^b- Vfe^ + V, 



0J3 



(9) 



(10) 



We label as Uk, ul and fk, fl the eigenvectors and eigenforms associated to the eigenvalues 
Afe, A^ = -Afc respectively, i.e., kuk = XkUk, Aul = -Xkul, fkA = Xkfk, /^A = -Xkfk, 
k — 1,2, 3, the * denoting complex conjugation. An explicit calculation leads to: 



U2 = S2 



/i = ti (Vb^ + V, iVb^ + V, 0, -i, 1, 0) , 
/2 = t2 (^-Vb^ + v, -iVb^ + V, 0, -i, 1, 0) , 



^^3 



= ss (0, 0, 0, 0, /3 = (0, 0, iV^, 0, 0, 1) , 
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where sj, tj G C, j = 1, 2, 3. We require that the eigenvectors and eigenforms are dual to each 
other [5,6, 13], namely, fjUk = f*ul = 6jk, fjul = f*Uk = 0, implying that si = ^,S2 = ^, 
S3 = The constants tj will be fixed later to simplify some commutation relationships. 
Finally, the eigenvectors and eigenforms satisfy the unit matrix decomposition 

3 3 
1 = XI ^^'^ ® ^'^ + ^fc ® ^fc) ^ A = XAfe(ufc®/fc-u*®/*) (11) 

k=l k=l 

® denoting tensor product. The A-expression in (fTT)) allows to decompose the Heisenberg 
trajectories as three oscillating modes of frequencies ujj [5,6]. Moreover, it will characterize as 
well the algebraic structure of the Hamiltonian. 



4 Algebraic structure of H 

We can define now three pairs of ladder operators of H, Lk = f^t], L\ = fkf], k = 1,2, 3, 
which obey the following commutation relations with H: 

[H,Lk] = ~tfl[tH,r,] = -oOkLk, [H,lI] = uj^lI (12) 

An explicit calculation leads to: 

Li = tl [V^T^iX - tY) + i{P, - iPy)] , 
L^ = tl[-y/WT^{X -tY)+t{P,-tPy)], L3 = tl{-tV^Z + P,). (13) 

By evaluating next the commutators between Li,L}-, the following non-null results are obtained: 

[Li,Ll] = 2|ti|2(^i-^2) = l, 
[L2,L\] = -2|t2pK -u;2) = -1, [^3,4] = 2|t3p^3 = 1, (14) 

where we have finally chosen tj G M"*" such that ti = t2 = 1/ s/^iuJi — UJ2), = 1/ \/2ujI to 
simplify at maximum equation ([141) . On the other hand, [L\, L^] = [Li, Lj] = 0, i, j = 1,2, 3. 
Now H is factorized in terms of Lk, L\ as follows [5, 6]: 

H = ujiL{Li - UJ2L2LI + UJ3LIL3 + {uji - UJ2 + uj3)/2. (15) 

Moreover, equations (I14I15I) allow to identify three independent oscillator modes for H, each 
one characterized by its number Nk, annihilation Bk and creation bI operator, in the way: 

Nk = BlBk, k = 1,2,3, (16) 
Bi=Li, B2 = L\, B^ = Ls, B\ = L\, bI = L2, bI = lI. (17) 

They obey the standard commutation relations: 

[Nk, Bk] = -Bk, [Nk, Bl] = Bl [B„ bI] = 5,k. j, k = 1,2, 3. (18) 
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Hence, one can construct a basis {|ni, n2, ns), rij = 0, 1, 2, ... ,j = 1, 2, 3} of common eigen- 
states of A^i, N2, N3, 

Nj\ni,n2,n3) = nj\ni,n2,n3), j = 1,2,3, (19) 
departing from an extremal state |0, 0, 0) which is annihilated by Bi, B2, B^: 

5,|0,0,0) = 0, J = 1,2,3. (20) 
If we assume that |0, 0, 0) is normalized, it turns out that [14]: 

\ni,n2,n3) = (^i! ^2! ng!)-!/^^!"! ^tnajo, Q, 0). (21) 
Moreover, Bj, Bj, j = 1, 2, 3 act onto \ni,n2,ns) in a standard way: 

Bi\ni,n2,n-i) = ^/ri7 1, ng, ng), BI \ni,n2,n3) = y/ni+ 1 \ni+ 1,^2,^3), 

and similar expressions for the action of B2, B\, B^, B\. Notice that 1^1,712,^3) is eigenstate 
of the Penning trap Hamiltonian with eigenvalue En^^n2,nz = 1^1(^1 + 1/2) — uj2{n2 + 1/2) + 
^^3('^3 + 1/2) = E{ni, n2, n^). In particular, the extremal state |0, 0, 0) has eigenvalue -£^0,0,0 = 
(cui —UJ2 + 1^3) /2, i.e., it is neither a ground nor a top state since its energy is "in the middle" of 
the spectrum of H. Following [15], it is seen that there is an intrinsic algebraic structure for our 
system, which is characterized by a linear relationship between the Penning trap Hamiltonian 
H and the three number operators A^^^: 

H = E{Ni, N2, N3) = uJiNi - UO2N2 + uj-iN^ + Eo,o,o- (22) 

As it happens for one-dimensional systems, in our three-dimensional example the detailed struc- 
ture is contained in the operator relation (|22l) . which is responsible of the specific spectrum and, 
consequently, of the lack of a ground or a top proper energy. On the other hand, the global struc- 
ture comes from the very existence of the three independent oscillator modes for H, each one 
characterized by the standard generators {Nj, Bj, j = 1, 2, 3. This global behavior al- 

lows to identify in a natural way the extremal state |0, 0, 0) which, although is neither a ground 
nor a top energy eigenstate, plays the same role as the ground state for the one-dimensional 
harmonic oscillator. 



5 Extremal state wave function 

The existence of the extremal state |0, 0, 0) is guaranteed by a theorem which is proven else- 
where [5]. It ensures that, if the operators 

Bj = iP-aj + R-f3j, B] = -ia]- P + ■ R, j = 1,2,3, (23) 

obey the commutation relations (fTSl ). then the system of partial differential equations {f\Bj\0, 0, 0) = 
0, j = 1, 2, 3, for the extremal state wave function 0o(r) = (^0; 0; 0) has a square integrable 
solution given by 

00 (^) = cexp (-^aijrirj) = cexp ( -^f'^af ] , (24) 
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where a = (aij) is a complex symmetric matrix satisfying 

aai=/3j-, j = 1,2,3. (25) 
According to (|23l) . through equations (|13I17I) we identify the vectors 

52 = P2 = ib' + v)'/'a,, (26) 



"3 



(0, 0, 1)^, 03 = i-2vy^^a3. 



V2{-2vy/^ 

Thus, a = diag [y/b"^ + f , v^P"H-t', V— 2f J , and from (|24|) we finally get the extremal state 
wave function we were looking for: 



Mrl = cexp I -^^^ix' + y^) - ] . (27) 



6 Penning trap coherent states 

Once the Penning trap Hamiltonian has been expressed appropriately in terms of annihilation 
and creation operators, we can develop a similar treatment as for the harmonic oscillator to build 
up the corresponding coherent states. 



6.1 Annihilation operator coherent states 

In the first place let us look for the annihilation operator coherent states (AOCS) as common 
eigenstates of Bi, B2, B-y. 

Bj\zi,Z2,Z3) = Zj\zi,Z2,Z3), j = 1, 2, 3. (28) 

Following a standard procedure, let us expand them in the basis {|ni, n2, n^)}: 

oo 

kl, 2^2, Z3) = ^ 

Cm ,712,713 

\ni,n2,n3). (29) 

ni,n2,n3=0 

By asking that (|28l ) is satisfied, three recurrence relationship for c„^_„2^„3 will be obtained, which 
in turn lead to the following expressions: 

Cni,n2,n3 = (^1 !) "^/^^^ C0,n2,n3 = ^-^^"^ Z2' Cn„0,n, = [m])-^/'^ z'^' Cn,,n2,0- (30) 

Hence, it is straightforward to show that 

Cnun„n, = {Ull Tla! ^3 !) "^/'^^ 4' Co,0,0, OD 

where co,o,o is to be found from the normalization condition. Thus, up to a global phase factor, 
the normalized AOCS become finally: 

|zi,^2,^3) = e-^^^^^^^^ f2 {n,\n2\n3\)-'/'z^'zrzr\n,,n2,n3). (32) 

rii,n2,n3=0 
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6.2 Displacement operator coherent states 

According to equation (O, for the j-th mode of the Penning trap Hamiltonian we have to take 
into account the corresponding displacement operator Dj{zj) = exp{zjBj — z*Bj). By using 
the BCH formula it turns out that: 

Dj{zj) = e-^e"^■^ie-"^*■^^ j = 1, 2, 3. (33) 
Now, the global displacement operator is given by: 

D(z) = D{zi, Z2, zs) = D,izi)D2iz2)Dsiz3), (34) 

where z denotes the complex variables zi, Z2, Z3 associated to the three modes. By employing 
now the second definition, we get the displacement operator coherent states (DOCS) |z) from 
applying -D(z) to the extremal state |0, 0, 0): 

I \ r,f Mn n n\ +1^2! +1^3! z^ Z2 z^ 712,^3; 

|z) = D(z)|0,0,0) = e 2 > . , , , • (35) 

n.i, 712,713=0 V i z d 

By comparing (l32l) and (l35l) ) we realize that the DOCS and the AOCS are the same. Moreover, 
since [z^B] - z*Bj, ZkBl - zlBk] = 0, j, k = 1,2, 3, we get 

D{z) =exp{ziBl + Z2BI + Z3BI - zlBi - z*B2 - zlB^i) =ex^[i{t-R - f-P)] 

= C{z)F{R) exp(-zf • P) = [C(z)]-i exp(-zf ■ P)F{R), (36) 



where we have used the BCH formula and equation (1231) to identify 

(6^ + t,)-i Re[zi - 2:2] \ / + v)^ lm[zi - Z2] 

-(6^ + v)-'^ lm[z, + Z2] \, S = (62 + Re[z, + Z2] | , (37) 
-{-v/2yi lm[z3] J V {-8v)-^Re[z,] 

C(z) = e-*^-^/2 ^ exp{z(Re[zi]Im[z2] + Re[z2]lm[zi] + Re[zs]lm[z3])} , 

F{R)=e'^-^= exp{i{b^+v)-i{lm[zi-Z2]X + Re[zi + Z2]Y)+i{-8v)-iRe[zs]Z}. 

Since the operator e~*^ '", Fj G R, performs a coordinate displacement in the way (r|e~*^'^ 
(r — r|, we finally get: 

,#>,(rl = (flz) = (r-lZ)(z)|0,0,0) = C(z)F(f){fle-"'-''|0,0,0), 

C(z)F(.>. I . - M£l^, , + M£l±^, , + (Z2) ' I . (38) 



with (/)o(r) given by (|27 
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7 Mean values of physical quantities 



Let us evaluate next the mean values (-Rj)z = (z|i?j|z), (Pj)z = (z|Pj|z), j = 1, 2, 3, and the 
corresponding mean square deviations in a given CS |z). To do that, we analyze first how the 
operators Rj, -R|, Pj, Pj are transformed under -D(z). By using equations (l36l) and (tST]) it is 
straightforward to show that: 

D\z)K;D{7) = {R,+V,Y, D\7)PfD{z) = {Pj + T.^)\ n = l,2,... (39) 

Therefore: 

(i?,)z = (i?,)o + r,-, = (i?2)^ + 2r,(i?,)o + r/, (Ai?,)^ = (Ai?,)^, (40) 

(P,)z = (P,)o + S„ {Pf), = {P^)o + 2S,(P,)o + S/, (AP,)^ = (AP,)^- (41) 

Notice that the mean square deviations of Rj and Pj are independent of Zi, Z2, but depend 
on {Rj)o, {Pj)o, {R'j)o, (-P/)o, j = 1, 2, 3, which need to be evaluated. The first six quantities 
can be obtained from the homogeneous equations (Pfc)o = ^(«fc)j (-P;)o + i/3k)j{Rj)o = 0, 
{bI)o = -iial)j{Pj)o + 0*k)j{Rj)o = 0, k = 1,2,3 (see (ES and use that Pfc|0,0,0) = 
(0,0,0|Pj = 0). By using (|26|) . the system to be solved becomes: 

V ¥T^i {X)o - t{Y)o) + t m)o - i{Py)o) = 0, 
-V^T^ ((X)o + i{Y)o) - I ((P.)o + i{Py)o) = 0, 

and the complex conjugate equations. Its solution is given by 

(P,)o = (P;)o = 0, j = 1,2,3. (42) 

In order to obtain {R^)o, {Pj)o, we calculate the mean values for the several products of pairs 
involving Bj, P^. From these thirty six equations just twenty one are linearly independent: 
{BjBk)o = 0, J = 1,2, 3, A; < j (six equations); (pJp1)o = 0, j = l,2,3,fc < j (six 
equations); {B\Bj)Q = 0, j,k = 1,2,3, (nine equations). By solving this linear system, the 
non-null mean values of the twenty one independent products of Pj and Pj are now: 

{X')o = {Y\ = [4(6^ + v)]-K {Z\ = {-SvrK 
{P')o = {P;)o = [{b' + v)/A]K (P|)o = {-v/2)K 
(XP,.)o = {YPy)o = {ZP.)o = t/2. 

The previous formulas imply that equations (|40I41I) become 

{AX)l = {AY)l = m' + vrK {AZ)l = {-8vrK 
{AP^)l = {APy)l = [{b' + v)/4, {AP,)l = {-v/2f^, 

and therefore 

(AX)z(AP,.)z = (Ar)z(AP,)z = (AZ)z(AP,)z = 1/2. 
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This means that our CS have minimum Heisenberg uncertainty relations. 

Finally, by using equations (I16I22I) we calculate the mean value of the Hamiltonian H in a 
given CS |z): 

(H)^ = - uj2\z2\'^ + cus^sP + ^o,o,o- (43) 

A similar calculation for {H'^)z can be done, leading to: 

{AH)l= (^b + Vb^y\zi\^+(b-V^T^y\z2\'~2v\z3\^. (44) 

Once again, the fact that H is not positive definite is clearly reflected in (l43l) . 

Along this work we have assumed that b = —4^ > 0. For 6 < 0, small differences concern- 

o 2mc ' 

ing the identification of the appropriate annihilation and creation operators arise. However, the 
extremal state and CS wave functions (poif^, 0z(^) as well as the corresponding mean values, 
will coincide with those previously calculated. In particular, the Heisenberg uncertainty relation 
will achieve once again its minimum value [14]. 

8 Concluding remarks 

In this paper it was introduced a technique to find the CS for a charged particle in an ideal 
Penning trap. We have shown that the coherent states, calculated through both definitions given 
by equations (|1I2|) . are the same. We introduced also a prescription to obtain the mean values of 
several physical observables in a given coherent state. We have found, finally, that the Penning 
trap coherent states (derived algebraically) obey also the third CS definition, i.e., they satisfy 
the minimum Heisenberg uncertainty relation. 

Let us remark that the method presented here is quite general, and it could be applied to other 
systems characterized by quadratic Hamiltonians. In order to implement systematically this 
treatment, we have to identify first the stability regions where the non-degenerate eigenvalues 
of A become purely imaginary, which ensures that the Heisenberg and classical trajectories are 
trapped. In the trap regime the Hamiltonian is decomposed in terms of independent harmonic 
oscillators, and thus our procedure can be straightforwardly applied. Note that generalizations 
of this kind have been elaborated elsewhere (see e.g. [16]). However, in our method it is direct 
to identify the global sign accompanying each individual oscillator involved in the Hamiltonian 
decomposition. As we saw in our Penning trap example, those signs determine the existence 
or not of a ground state for the system, a fact which is not well known in the literature. More- 
over, they become fundamental for the determination of the intrinsic algebraic structure of the 
involved Hamiltonian (compare equation (l22l)). Observe that some of these properties were 
found previously for operators imitating the Hamiltonian in non-inertial reference frames [5,6]. 
By means of this example we have shown that such a property arises as well for Hamiltonians 
in inertial frames of reference. 
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